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ABSTRACT 



Context. The Cygnus region hosts a giant molecular-cloud complex that actively forms massive stars. Interactions of cosmic rays with interstellar 
gas and radiation fields make it shine at y-ray energies. Several y-ray pulsars and other energetic sources are seen in this direction. 
Aims. In this paper we analyze the y-ray emission measured by the Fermi Large Area Telescope in the energy range from 100 MeV to 100 GeV 
in order to probe the gas and cosmic-ray content on the scale of the whole Cygnus complex. The y-ray emission on the scale of the central massive 
stellar clusters and from individual sources is addressed elsewhere. 
^-H Mettiods. The signal from bright pulsars is greatly reduced by selecting photons in their off'-pulse phase intervals. We compare the diff'use y-ray 
emission with interstellar gas maps derived from radio/mm- wave lines and visual extinction data. A general model of the region, including other 
pulsars and y-ray sources, is sought. 

r* ' Results. The integral H i emissivity above 100 MeV averaged over the whole Cygnus complex amounts to [2.06 + 0.11 (stat.) (syst.)] x 10~^^ 
, ^ photons s~^ sr~^ H-atom~\ where the systematic error is dominated by the uncertainty on the Hi opacity to calculate its column densities. The 
integral emissivity and its spectral energy distribution are both consistent within the systematics with LAT measurements in the interstellar space 
^ near the solar system. The average Xco = A^(H2)/Wco ratio is found to be [1.68 + 0.05 (stat.) +q^q(Hi opacity)] x 10^° molecules cm~^ (K 
km s~^)~\ consistent with other LAT measurements in the Local Arm. We detect significant y-ray emission from dark neutral gas for a mass 
' corresponding to ~ 40% of what is traced by CO. The total interstellar mass in the Cygnus complex inferred from its y-ray emission amounts to 
8 +^ X 10^ Mo at a distance of 1.4 kpc. 

Conclusions. Despite the conspicuous star formation activity and high masses of the interstellar clouds, the cosmic-ray population in the Cygnus 
complex averaged over a few hundred parsecs is similar to that of the local interstellar space. 

Key words. ISM: abundances - ISM: clouds - cosmic rays - Gamma rays: ISM 



1. Introduction 

Regions with conspicuous star formation activity are of great in- 
terest for understanding the life cycle of interstellar matter and 
the properties of cosmic rays (CRs) in the Galaxy. Interstellar 
y-ray emission produced by CR interactions with the interstel- 



lar gas via nucleon-nucleon inelastic collisions and electron 
Bremsstrahlung can be used to probe their CR and gas content. 

High-energy y-ray observations have entered a new era since 
the launch of the Fermi Gamma-ray Space Telescope in 2008. 
The Fermi Large Area Telescope (LAT: lAtwood et al l[2QQ9) has 
already measured strong y-ray emission toward the 30 Doradus 
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starburst region in the Large Magellanic Cloud (lAbdo et alJ 
l2010dh . and it also pointed out a global correlation between the 
y-r ay luminosity and s tar-formation rate in a few normal galax- 
ies (lAbdo et al.ll2010d) . 

A primary observational target for Fermi in our Galaxy is 
the Cyg nus X star-forming region, owing to it s proximity (~ 
1.4 kpc; lHanson|[2003l: iNegueruela et alJ 12008') and the avail- 
ability of numerous multiwavelength observ ations. Named af - 
ter the strong emission at X-ray wavelengths (ICash et al.ll 19801) , 
Cygnus X is located around the Galactic longitude / = 80°, tan- 
gent to the Lo cal Spur. It cont a ins nu n ierous H ii regions and OB 
associations (lUvaniker et al.l 1200 ll: |lE Duigou & Knodlsedej 
l2002h . It has long been debated whether it represents a coherent 
complex or the alignment of different structure s along the line 
of sigh t. Re cent h igh-resolution observations by lSchneider et al.l 
(l2006h and lRov et a l. ( 2011) have pointed out that most of the 
molecular clouds in the Cygnus X region are connected and 
partly show evidence of interactions with the massive stellar 
cluster Cygnus 0B2 and other OB associations in the region. 
Foreground molecular clouds from the Great Cygnus Rift, at 
0.6-0.8 kpc, contribute little to the high mass seen in interac- 
tion with the Cygnus X region itself, at 1.4 kpc. Therefore, the 
molecular cloud complex appears as one of the most massive in 
the Galaxy. Atomic gas seen in these directions is probably more 
widespread along the line of sight. 

lAbdo et al.l (l2007l l2008l) analyzed Milagro measurements at 
energies > 10 TeV and reports an excess of diffuse y-ray emis- 
sion with respect to predictions based on CR spectra equivalent 
to those near the Earth. They attribute the excess to the possible 
presence of freshly accelerated particles. 

The escape of CRs from their sources and the early propa- 
gation in the surrounding medium have so far been poorly con- 
strained by observations. In particular, particles accelerated in 
regions of massive star-formation are likely to be significantly 
influenced by the turbulent environment. It is therefore interest- 
ing to investigate how the CR populations on the scale of the 
massive stellar clusters and on the larger scale of the parent in- 
terstellar complex compare with each other and with the average 
CR population of the Local Spur. 

This paper reports our analysis of the y-ray emission mea- 
sured by the LAT in the energy range between 100 MeV and 
100 GeV across the entire Cygnus region. We focus on the large- 
scale properties of the interstellar emission to probe the CR pop- 
ulation and to complement gas and dust observations at other 
wavelengths to constrain the amount of gas in different phases 
over the whole Cygnus complex. We also build an improved in- 
terstellar background framework for the study of individual y- 
ray sources that will be treated in companion papers. We discuss 
interstellar emission in t he star-forming re gion of Cygnus X in a 
dedicated paper (lAckermann et al.|[201 lah . 



2. Data 

In this section we describe the data used in the paper. First of all, 
we describe the y-ray data sample (§ 12. lb , including the selection 
criteria (§ 12.1.11) and the procedure to mitigate the signals from 
the brightest pulsars (§ 12.1.21) . Then, we present an overview of 
the other data used to trace the distribution of the interstellar mat- 
ter (§ 12.21) , including radio and mm- wave lines (§ 12.2.11) to trace 
neutral gas, visual extinction (§ 12.2.2b to trace dark neutral gas, 
and the free-free emission intensities obtained from microwave 
observations (§ 12.2.31) to trace the ionized gas. 



2.1. Gamma-ray data 

2.1 .1 . Observations and data selection 

The LAT is a pair-tracking telescope detecting photons from 
20 MeV to inore tha n 300 GeV. The instrument is descri bed in 
Atwood et al.l (I2009h and its on-orbit calibration in Abdo et all 
(I2009ah . The LAT operates most of the time in continuous sky- 
survey mode. We accumulated data for our region of interest 
from August 5, 2008 (METQ 239587201) to August 5, 2010 
(MET 302659202). 

We selected data according to the tightest available back- 
ground reje ction criteria, corres ponding to the Pass 6 Dataclean 
event class (I Abdo et al.ll2010elPl . In order to limit the contami- 
nation from the Earth's atmospheric y-ray emission, we selected 
events with measured arrival directions within 100° of the local 
zenith and within 65° of the instrument boresight, taken during 
periods when the LAT rocking angle was less than 52°. 

The angular resolution of the LAT strongly d epends on the 
photo n energy, improving as the energy increases (lAtwood et al.l 
[2009h . Confusion at low energies is a problem since we aim 
to spatially separate the different components in the crowded 
Cygnus X region. We therefore accepted below 1 GeV only pho- 
tons that produced electron-positron pairs in the thin converter 
planes of the front s ection of the tracker , which provides a higher 
angular resolution (I Atwood et al.ll2009l) . Above 1 GeV, we kept 
all events which converted either in the front or back section of 
the tracker. 

We analysed data at Galactic longitudes 72° < / < 88° and 
latitudes -15° < b < -1-15°. The longitude window contains 
the interstellar complexes associated with Cygnus X; the lati- 
tude window is large enough to allow a reliable separation of 
the large-scale emission from atomic gas, isotropic background 
and inverse-Compton (IC) scattering of low-energy radiation 
fields by CR electrons. We analysed the data in the 100 MeV- 
100 GeV energ y band. Below 100 MeV the instrumental system- 
atics are large (iRando et and the angular resolution is 

poor, whereas above 100 GeV we are limited by the low photon 
statistics. 

2.1 .2. Removal of bright pulsars 

Three bright pulsars dominate the y-ray emission f rom the re- 
gion b elow a few GeV: the radio pulsar J2021 -h365 1 (lAbdo et al.l 
l2009el) and t he two LAT-discoy ered pulsars J202 1-^4026 and 
J2032-F4127 (lAbdo et al.ll2009bh . To increase the sensitivity to 
faint sources and to the spatial structure of the diffuse emis- 
sion, we reduced their contribution by excluding the periodic 
time intervals when their pulsed emission is the most intense. 
Removing the intense pulsed flux helps to reduce the impact of 
any incorrect modeling of such bright sources on the results. 

To assign pulse phases for each of the three pulsars , we pro- 
duced timing models using T empo2 (iHobbs et al.ll2006h accord- 
ing to the method described in lRay et al.l (1201 lljC Figurefllshows 
the three light curves and the phase intervals with bright pulsed 
emission. The phase boundaries are reported in Tabled together 
with the fraction of time in the off-pulse interval suitable for our 
study. There is a considerable level of off-pulse emission toward 

^ Fermi Mission Elapsed Time, i.e. seconds since 2001 January 1 at 
00:00:00 UTC. 

^ Performance figures for the Dataclean event selection are given in 
the reference. 

^ For the three pulsars, the RMS of the timing residuals is below 1.1% 
of their rotational period. 
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Fig. 1. Light curves and ofF-pulse phase intervals for the three 
bright pulsars. The light curves are constructed for illustration 
purposes with photons recorded in a circular region of radius 
0.5° around the pulsar positions and energies > 200 MeV. 

Table 1. Off-pulse phase intervals and time fractions 



of the three bright pulsars. 


PSR 


phase interval 


time fraction (%) 


J2021+3651 


0-0.12, 0.22-0.59, 0.7-1 


79 


J2021+4026 


0.18-0.5, 0.62-0.94 


64 


J2032+4127 


0-0.14, 0.2-0.59, 0.79-1 


74 



PSR J202 1+4026 that cannot be removed (lAbdo et al.ll2010bh : 
however, given the brightness of the source, the removal of the 
on-pulse interval is useful for our aims. 

A total count map in the off-pulse phase intervals of the three 
bright pulsars is provided for illustration in Fig.O To remove the 
pulsar signal without excessively sacrificing the photon statistics 
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Fig. 2. Total count map in the energy range 100 MeV- 100 GeV, 
binned over a 0.125° x 0.125° grid in Galactic coordinates in 
Cartesian projection. Data were selected according to the crite- 
ria described in the text (§ 12.1.11) and in the off-pulse phase in- 
tervals of the three bright pulsars (§ 12.1.21) , whose positions are 
marked by diamonds. Counts are saturated between and 50, 
and smoothed for display with a Gaussian kernel of cr = 0.25°. 



in other directions, we restricted the timing selection to a circular 
region around the pulsar position, namely to pixels in our angular 
grid (described later in § 13.21) the centroids of which lie within 
the energy-dependent radius 



^cut(^) = 2 • 



0.8° (- 



GeV/ 



-0.8 



e 0.07^ 



(1) 



where the symbol © indicates addition in quadrature. This is an 
approximate representation of the LAT 95% containment angle 
as a function of energy. We note that the accurate parametriza- 
tion of the LAT point spread function (PSF) depends on energy, 
pair-conversion point in the tracker and, to a lesser extent, on 
the incidence angle. The PSF is best represented by the LAT in- 
strument response functions (IRFs), which are used later for the 
likelihood analysis. The above acceptance-averaged approxima- 



tion for the containment angle is only useful for calculating the 
radius rcut, and we verified that the results are insensitive to rea- 
sonable variations in this parameter. 

To take the cut on pulsar phases into account, for each di- 
rection in the sky and energy the exposure (see again § 13.21) was 
multiplied by the remaining livetime fraction. The remainder of 
the pulsar emission was included in the model using 

- a point source to represent emission in the off'-pulse interval; 

- a second point source, for which the number of expected 
counts is set to zero at r < rcut(E) from the pulsar position, 
to represent on-pulse y-rays spilling over at r > rcut(E) in the 
tails of the PSF. 

The two sources have free independent fluxes in each energy bin 
of the analysis. This is particularly important to account for the 
diff'erent spectra of the on-pulse and ofl'-pulse y-ray emission and 
also to compensate for any mismatch between the tails of the 
model PSF and the emission at large angles from the brightest 
sources in the region. 

Since th e three pulsars have exponential spectral cutofl's near 
2-3 GeV dAbdo et al.ll2010bl) the phase selection was not ap- 
plied above 10 GeV where the level of pulsed emission is low 
and each pulsar was accounted for by a single point source. On 
the other hand, given the abundant statistics but poor angular 
resolution at low energies (more than half of the region of in- 
terest would be subject to on-pulse event removal), we selected 
ofl'-pulse photons for the whole region below 316 MeV0. In this 
case no "spill-over" source was necessary. 

2.2. Ancillary data 

2.2.1 . Radio/mm-wave lines: neutral gas 

Neutral atomic hydrogen. Hi, was traced thanks to its 21- 
cm line. Where available0 w e used data froni the Canadian 
Galactic Plane Survey rCGPS: Pl^lor et al.|[2003b rebinned onto 
the 0.125° X 0.125° grid used for the other maps. Elsewhere, we 
used data from the Leiden/ Argentine/Bonn (LAB: [Kalberla et aD 
120051) survey, with a coarser binning of 0.5°. We checked the 
consistency of the two survey calibrations in the overlap region. 

Molecular hydrogen cannot be observed directly in its most 
abundant cold phase. The velocity-integrated brightness temper- 
ature of the ^^CO 2.6-mm line, Wqo, is often assumed to lin- 
early scale with the N(}l2) column densit y. We used CO data 
from the composite survey by Dame et'aP (l200ll) , filtered with 
the moment-masking technique ( Damdl2011l) ' in order to reduce 
the noise while preserving the faint cloud edges and keeping the 
resolution of the original data. 

The Doppler shift of radio/mm-wave lines can be used to 
kinematically separate the Cygnus complex from two faint seg- 
ments of the Perseus and outer spiral arms seen beyond Cygnus 
in the same directio n. We applied the k inematic separation pro- 
cedure illustrated by lAbdo et al.l (l2010j) . starting from a prelim- 
inary boundary located at a Galactocentric radiu^T? = 9.4 kpc 
and then adapting the separation to the cloud structures and cor- 
recting for the spill-over due to the broad velocity dispersion of 
Hi lines. The separation into two regions is accurate enough to 



model the in terstellar y-ray eni i ssion in Cygnus since lAbdo et al.l 
(1201 Oft and lAckerma nn et al.l (l2011bb did not find any signifi- 
cant gradients of the gas y-ray emissivities in the outer region of 
the Milky Way. We applied the kinematic separation procedure 
to prepare maps of the column densities of atomic gas, N(tl l), 
and of Wco- The maps are shown in Fig. [3] for H i and Fig. |4] for 
CO. All the gas maps mentioned in the paper have > 10° borders 
around the analysis region used to properly convolve the model 
with the LAT PSF. 

Substantial uncertainties in the determination of A/^(H l) arise 
from the choice of spin temperature for the optical depth correc- 
tion. We adopted a uniform Ts = 250 K as baseline case, which 
is the average spin temperature that best reproduces the blend- 
ing of cold and warm atomic gas according to observations of 
em ission-absorption H I pairs in the region covered by the CGPS 
(Di ckev et al.ll2009h . Other values 100 K < < oo will be con- 
sidered later to evaluate the related systematic uncertainties af- 
fecting the results of our analysis. 

2.2.2. Visual extinction: dark neutral gas 

Multiwavelength observations indicate that the combination of 
the Hi and CO lines does not properly trace the total col- 
umn densities of the neutral inters t ellar medium (ISM) (e.g . 
"Magnani et al. 2003; Gre nier et al] 120051: lAbdo et al.l l2010fl 
Langeret al. 2010: Ade etaP l201ll) . Since the work by 
Grenier et al.l (l2005l) . dust tracers have been used in y-ray anal- 



^ See § l3.2| for the definition of the energy grid used in the analysis. 

^ The CGPS coverage is almost complete at -3.5° < b < +5.5° for 
this longitude range. 

^ Based on the assumption of a flat rotation curve with solar radius 
RQ - 8.5 kpc and Galactic rotation velocity at the solar circle VO = 
220 km s-^ 



yses to complement the H i and CO lines, under the assumption 
that dust grains are well mixed with gas in the warm and cold 
phases of the ISM and therefore provid e an estimate of total gas 
column densities. iGrenier et aP (l2005l) a ndlAbdoetan (1201 Of|) 
adopted the E(B - V) color excess map by lSchlegel et a l. (1998) 
as a tracer of the total column densities, and used the E(B - V) 
residuals -i.e. E(B - V) minus the best-fit linear combination of 
A/^(H l) and Wqo maps- as a tracer of the dark-gas column densi- 
ties in nearby clouds of the Gould Belt. 

The use of the E(B - V) map is problematic in the Cygnus X 
region for two reasons: 

- numerous infrared point sources contaminate the map; 

- the temperature correction used by Schlegel et'aP (Il998l) to 
derive the dust column-density map from IRASI\SSA mea- 
surements is highly uncertain in regions of massive star- 
formation because of the enhanced radiation fields. 

We have therefore adopted the visual extinction Ay as 
derived from th e reddening of near -infrared sources in the 
2 MASS catalog (ISkmtskie et al.ll2006l). The Ay maps produced 
bv lRowles & FroebrichI (l2009h and iFroebrich & RowlesI (l2010h 
were used for Ay < 5 mag. They exhibit saturation eff'ects 
at higher extinction values, so we complemented them with 
an Ay map obtained from 2MASS d a ta usin g the code and 
method developed by ISchneider et al.l (1201 ll). The latter use 
the Besangon stel lar population model (iRobin & Crezd Il986l: 
iRobin et al.l I2OO3I) to filter out the contribution from the fore- 
ground bluest starfl. The second Ay map was built in a 12° 
region centered on (/, Z?) = (80°, 0°), and, compared with the 
first set of maps, it presented an off'set of ~ 0.46 mag at low 
extinction. We cons t ructed the final Ay map from the direct 
iRowles & FroebrichI (|2009') data below 5 mag and from the sec- 
ond map, off'set by 0.46 mag, at higher extinction. 



^ To do so, a distance from the observer needs to be assumed for the 
clouds under consideration. We verified that variations of a few hundred 
parsecs do not significantly change the results presented in the paper. 
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Fig. 3. Maps of NiYL l) column densities in the Cygnus complex in the Local Spur (left) and in the outer Galaxy (right), under the 
assumption of a uniform spin temperature of 250 K. The color scales with 7V(Hl) in units of 10^^ atoms cm"^. The maps were 
smoothed with a Gaussian kernel of cr = 0.25° for display. 



The Ay map was binned onto the same 0.125° x 0.125° grid 
in Cartesian projection as the other maps. The Ay map was fit- 
ted with a linear combination of the A/^(H l) and Wqo maps pre- 
viously described. The input Ay map minus the best-fit linear 
combination of the A/^(H l) and Wqo maps yielded the Ay excess 
map, Ay^exc, which will be used to trace the dark neutral gas. 
Only residuals corresponding to input Ay > 0.3 mag were kept 
to limit the noise off' the plane. The Ay excess map is shown in 
Fig.E 



2.2.3. Microwave emission: ionized gas 



Away from H ll regions around massive stars and stellar clus- 
ters, the ionized gas constitutes a layer of characteristic height 
> 1 kpc over the Galactic plane with lit tle mass compared to 
the neutral phases (ICordes & Lazioll2QQ2l) . Therefore, it has of- 
ten been neglected in previous y-ray studies. However, we find 
in the Cygnus X region many conspicuous H 1 1 regions excited 
by the intense radiat i on fields of the nu merous massive stars 
(lUvaniker et al.ll200ll:IParadini et al.ll2003h . 



Ionized gas masses ca n be traced by free-free emi ssion fol- 
lowing the prescription by lSodroski et al.l (fl 98911 199?!) to derive 
the N(ll II) column densities: 

^(Hii) = 1.2xl0-cm-(i^)"t-^ry-^f '-^,(2) 
\1K/ \lcm-3/ VlGHz/ Uysr"! 

where /ff is the free-free emission intensity at the frequency v, 
is the electron temperature, and ^eff the eff'ective electron number 
density. We adopted a free-free emission map deriv ed from the 
seven- year WMAP data in the Q band (40 GHz) bv lGold et aP 
(1201 Ih using the maximum entropy method from the prior tem- 
plate given by the extinction-corrected Ha map by Hnkbeinej 
(2003). It was rebinned onto the 0.125° x 0.125° grid used for 
the other maps, as shown in Fig.[6j 

3. Analysis 

3.1. Analysis model 
3.1 .1 . Diffuse emission 

Since the bulk of Galactic CRs in the relevant energy ranges 
are expected to be smoothly distributed on scales exceeding 
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Fig. 4. Maps of Wqo intensities in the Cygnus complex in the Local Spur (left) and in the outer Galaxy (right). The color scales with 
Wco in units of K km s"^ above 1.5 K km s"^ The maps were smoothed with a Gaussian kernel of cr = 0.25° for display. 



the typical dimensions of interstellar clouds and to penetrate all 
the phases of the ISM uniformly, the y-ray emission produced 
by CR-gas interactions can be modeled to first order as a linear 
combination of the gas column densities summed for the diff'er- 
ent phases and diff'erent regions along the line of sight. 

Ionized gas, with a total mass over the region of 0.4 x 
10^ (^eff/1 cm"^)-^ Mo for Te = 10"^ K, represents less than 4% 
of the total atomic mas s present in the Cygn us complex (assum- 
ing riQff = 2-10 cm'^. ISodroski et aO 19971) . The corresponding 
column densities are highest in the massive star-forming region 
of Cygnus X where we detected a bright and hard extended y- 
ray source powered by freshly-ac celerated particles, that w ill be 
called hereinafter "the cocoon" (lAckermann et al.|[201 lah . The 
free-free emission map was significantly detected in addition to 
the other interstellar components, but only at the expense of an 
unusually hard spectrum. To model the entire region, we intro- 
duced an extended source to account for the cocoon, as described 
in l3.1.2[ The latter was found to provide the best fit to the LAT 
data, yielding a higher maximum-likelihood value than the free- 
free emission template. Since the cocoon source overlaps most 
of the ionized clouds, it absorbs their contribution to the y-ray 
emission, so ionized gas was not included in the baseline model 
through the free-free emission template. The Cygnus X region is 



treated in detail in a companion paper (lAckermann et al.ll201 lah . 
The results presented in this paper were checked against the in- 
clusion of the free-free emission template in the model. 

The interstellar IC emission is produced by interactions of 
CR electrons and positrons with the low-energy interstellar ra- 
diation field (ISRF). To account for large-scale IC emission 
from the Milky Way we adopted a template calculated using the 
GALPROP CR pro pagation cod^ ( Strong & Moskalenko 1998; 
IStrong etal] 120071) . run 54_87Xexph7S. The IC emission was 
calculated on the basis of a CR electr on spectrum consis tent with 
recent measurements at the Earth (lAbdo et aDl2009cl) and the 
new calculation of the Galactic ISRF bv lPorter et al .1(1200 8h . 

Local radiation fields could leave unmodeled structures in IC 
emissio n, notably in the niassive star-forming region of Cygnus 
X (e.g. lOrlando & Strongll2007l) . In the companion paper we 
show that an upper bound to the IC emission from the stellar 
and interstellar low-energy radiation fields upscattered by CR 
electrons with the local spectrum is two orders of magnitude 
fainter than the cocoon emission, which in turn is fainter than the 
emission from the neutral gas (Fig.O. The CR electron sources 
within Cygnus X could further enhance the IC y-ray yield. Any 



http://galprop. stanford.edu 
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Fig. 5. Av excess map (magnitudes) obtained from the optical 
extinction Ay estimated from 2MASS data minus the best-fit 
linear combination of the A/^(H l) and Wqo maps shown in Figs. [3] 
and m respectively. The map was smoothed using a Gaussian 
kernel of cr = 0.25° for display. 



enhanced IC contribution from the inner region is accounted for 
in this analysis by the extended cocoon source, and it should not 
bias the determination of the gas emissivities we aim to study 
here. 

The diff'use emission model is completed by the isotropic 
background which combines the residual backgrounds from 
misclassified CR interactions in the LAT and the isotropic, 
presumably extrag alactic, y-ray emission (studied in detail in 
lAbdo et al.i[2010el) . 



3.1.2. Sources 

We included in the mod el the identified so urces in the region 
of interest: Cygnus X-3 (lAbdo et al.ll2009l) . PSR J1957+5033 
dSaz P arkinson et al.ll201Ql) and PSR J2030+3641 (ICamilo et al.l 
I2OIII) . in addition to the three bright pulsars as discussed in 
§[2T2l 
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Fig. 6. Free-free emission intensities from WMAP data. The 
color scales with brightness temperature in mK. The map was 
smoothed using a Gaussian kernel of cr = 0.25° for display. 



We also iterativ ely included significant IFGL point sources 
(lAbdo etal.1 l2010ah either associated with active galactic nu- 
clei (AGN) or characterized by variability or both. The 
sources were added with decreasing brightness: J2116.1-h3338, 
J2001.1-h4351, J2027.6-F3335, J2115.5-h2937, J2015.7-h3708, 
J2029.2-h4924, J20 12.2-^4629, and J2128.0-h3623. The iterative 
procedure is useful for stabilizing the likelihood fitting proce- 
dure and assessing the significance of sources added at each step. 

We detected extended y-ray emission above the global 
interstellar emission model discussed here associated with 
the supernova remnants known as the Cygnus Loop (G74.0- 
8.5 , e.g. Suneta l. 2006 ) and y Cygni (G78.2-h2.1, e.g. 
iLadouceur & Pineaultl l2008h and with the inner 100 pc of the 
Cygnus X complex. They are discussed in detail elsewhere. We 
briefly summarize here how these extended sources are modeled. 
For each of them we have tested difl'erent models and therefore 
verified that their presence does not bias the results concerning 
the properties of large-scale interstellar emission presented in the 
paper. 



The Cygnus Loop was modeled using a ring centered at 
(l,b) = (74.1°, -8.5°) and with inner/outer radii of 0.7° and 
1.6°, respectively, which best reprodu ces y-ray emission from 
the Cygnus Loop (iKatagiri et aLlboill) . 

We included two sources in the region of y Cygni in addition 
to PSR J2021+4026: 

- a uniform disk centered at (/, b) = (78.2°, +2. 1°) and a radius 
of 0.5° (G78.2+2.1; Green'!2009|); 

- a 2D Gaussian correspondin g to the moderately e xtended 
TeV soured] VER 2019+407 (IWeinstein et aDl2009l) . 

We detected extended y-ray emission toward the inner ~ 
100 pc of Cygnus X, which is effectively treated here as a source 
named "the cocoon". We discuss the nature of the cocoon and its 
relation with CR acc eleration in the massive st ar-forming region 
in a dedicated paper (lAckermann et al.|[201 lah . where we deter- 
mine a Gaussian centered at (l,b) = (79.6° ± 0.3°, 1.4° ± 0.4°) 
with a (T = 2.0° ± 0.2° width to be the model providing the best 
fit to the LAT data. As noted above, the cocoon source eff'ec- 
tively accounts for the contribution from ionized gas and from 
enhanced IC emission in the Cygnus X massive star-forming re- 
gion, as well as for locally-accelerated CRs. The spatial distribu- 
tion of the sources included in the analysis model is summarized 
in Fig. [71 

3.1 .3. Summary of the analysis model 

To summarize, the y-ray intensities / (photons cm"^ s"^ sr"^) are 
modeled in each energy bin by 

2 

1(1, b) = ^[^HM-A^(Hi)(/,Z7), + ^co,.- Wco(/,/^).] + 

1=1 

+ ^Av • Av,exc(/, b) -h IC(/, b) -h /iso + 

+ ^^//,^). (3) 

J 

The sum over i represents the combination of the two regions: 
1) Cygnus complex and 2) outer Galaxy. The free parameters 
of the diffuse emission model are the emissivities per hydrogen 
atom, qui,i (s"^ sr"^), the emissivities per unit of Wqo intensity, 
qco,i (cm"^ s"^ sr"^ (K km s"^)"^), the emissivity per Ay excess 
unit qAy (cm"^ s"^ sr"^ mag"^), and the isotropic intensity /iso 
(cm"^ s"^ sr"^). IC stands for the IC emission model described 
above. The sum over j represents the combination of the sources, 
either point-like or extended, as described in § 13.1.21 including 
a free parameter (flux normalization) independently for each of 
them. 
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Fig. 7. Sources included in the analysis model. Diamonds mark 
the positions of the three bright pulsars that were dimmed 
by phase selection (§ 12.1.21) . X points mark the positions of 
other identified sources. Crosses correspond to IFGL sources 
either associated to AGN or variable, or both, confirmed by our 
analysis. The blue circles correspond to the rims of the tem- 
plates adopted to model SNRs, the Cygnus Loop (G74.0-8.5) 
and y Cygni (G78.2-h2.1). The red circle marks the centroid 
of VER 2019+407 (whose extension is not appreciable in this 
large-scale view). The magenta circle represents the Icr contour 
of the Gaussian used to model the cocoon. 



3.2. Analysis method 

The model was fit to LAT data by using a binned maximum like- 
lihood with Poisson statistic^ independently over several en- 
ergy bins. We used a 0.125° x 0.125° binning in Cartesian pro- 
jection, comparable to the LAT angular resolution at the highest 
energies. We considered three energy bands: low (100 MeV-1 
GeV), mid (1 GeV-10 GeV) and high energies (10 GeV-100 
GeV). The low and mid-energy bands were divided further into 
four logarithmic- spaced energy bins, the higher-energy band in 

9 The 2D Gaussian fitted to TeV data by 'Weinstein et al. ( 2009) pro- 
vides a better fit to LAT data with respect to the coincident point source 
IFGL J2020+4049. 

As implemented in the standard LAT analysis tools 09-18-05. 



two because of the limited statistic£3- The analysis was based 
on the post-launch IRFs of the P6_V3 series, which consider ef- 
ficiency losses due to pile-up and accidental coincidence effects 
in the detector (iRando et al.ll2009l) . 

To perform the convolution with the LAT PSF, a power-law 
spectrum with index 2. 1 was assumed for the gas maps and other 
sources modeled by geometrical templates (the results do not 
significantly depend on this value). For all other sources we 
used power-law spectra with the spectral index reported in the 
IF GL Catalog. For the pulsars included in the LAT pulsar cata- 
log (lAbdo et al.l2010bh , we used the spectral functions described 
therein. 

The bounds of the energy bins are reported in Table [2 
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Fig. 8. Residuals (data-model). Left: low energies (100 MeV-1 GeV); center: mid-energies (1 GeV-10 GeV); right: high energies 
(10 GeV- 100 GeV). Units are approximate standard deviations (square root of model counts) saturated between ±3 and smoothed 
for display with a Gaussian kernel of cr = 0.25°. In each panel the blue circle in the top left corner represents the effective LAT 
PSF 68% containment circle for the event selection used in the analysis (averaged over the corresponding energy range assuming 
a power-law spectrum with index 2.1). Diamonds mark the positions of bright pulsars for which phase selection was applied, as in 
Fig. [21 Crosses mark the positions of unassociated IFGL sources coincident with positive residuals; the X point mark the position 
of 2FGL J2018.0 -h 3626, coincident with a hot spot in the residual map 1 GeV-10 GeV. 

Table 2. Best-fit parameters describing emission from interstellar gas (Eq.O under the assumption of a uniform H i spin temperature 
Ts = 250 K. 



energy bin"^ 


b 

^Hl,l 


^CO, 1 


b 

^Hl,2 


^C0,2 


d 

^Av 


0.1-0.178 


7.9+1.1 


3.3+1.0 


8.5+1.4 


0.00+0.06 


10+30 


0.178-0.316 


5.9+0.3 


2.2+0.3 


4.7+0.4 


0.000+0.002 


19+5 


0.316-0.562 


3.27+0.11 


1.16+0.08 


3.23+0.14 


0.00+0.06 


11.0+1.8 


0.562-1 


1.95+0.06 


0.59+0.04 


1.72+0.10 


0.5+0.3 


6.2+0.7 


1-1.78 


0.98+0.03 


0.328+0.016 


0.74+0.04 


0.12+0.13 


2.6+0.3 


1.78-3.16 


0.389+0.016 


0.141+0.008 


0.36+0.02 


0.02+0.06 


0.86+0.15 


3.16-5.62 


0.151+0.005 


0.044+0.004 


0.113+0.013 


0.02+0.03 


0.39+0.08 


5.62-10 


0.050+0.003 


0.016+0.002 


0.046+0.006 


0.000+0.005 


0.15+0.04 


10-31.6 


0.0085+0.0015 


0.0059+0.0010 


0.021+0.003 


0.002+0.008 


0.043+0.019 


31.6-100 


0.0024+0.0007 


0.0016+0.0004 


0.0007+0.0014 


0.002+0.003 


0.002+0.007 



Notes. Subscripts refer to the two regions separated in analysis: 1) the Cygnus complex in the Local Spur, 2) the outer Galaxy. Some parameters 
are poorly determined but they are reported for completeness. 

(a) (b) iQ-n g-i (c) ;^o-6 cm-2 s-1 sr-1 (Kkm s-i)-i 10"^ cm'^ s'^ sr'^ mag-^ 



4. Results and discussion 

4.1. Summary of the results and uncertainties 

The y-ray residuals corresponding to the best-fit model are 
shown in Fig. [8l They indicate that the model satisfactorily re- 
produces the morphology of the y-ray emission on larger angular 
scales than the LAT PSF in all the energy bands. Localized posi- 
tive residuals are still present. Some of them coincide with unas- 
sociated IFGL sources and others are associated with sources of 



the 2FGL catalo£3, notably 2FGL J2018.0 -f-3626also coinci- 
dent with the TeV source MGRO J2019-F27 (lAbdo et al.ll2007b . 
We verified that including sources accounting for those residu- 
als in the analysis model would not significantly aff'ect the deter- 
mination of y-ray emissivities associated with the diff'erent gas 
components summarized in Table [2j 



The preparation of the 2FGL catalog ran in parallel with the anal- 
ysis reported in this paper. The source list is now available from 
jhttp://ferm i.gsfc.nas a.goY/ssc/data/access/lat/2yr_catalog/^ 
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Fig. 11. H I emissivity spectrum in the Cygnus complex. Points: 
the best-fit estimate for the spin temperature Ts = 250 K. 
Hatched rectangles: systematic uncertainties taking H I opacity 
and y-ray selection efficienc y into account. Lin e: model of the 
local interstellar spectrum by 'Abdo et al.l l2Q09dl (with a nuclear 
enhancement factor of 1.84. iMori.2QQ9i) . 

and has a lower characteristic height above the Galactic plane, 
and can therefore be reliably determined in the likelihood fit. We 
verified that completely neglecting the large-scale IC emission 
leads to negligible variations in the emissivities of CO-bright 
gas and Ay excesses, and to variations lower than or comparable 
to statistical uncertainties for the emissivities of the atomic hy- 
drogen, which is less structured than the other ISM components 
and has a larger characteristic height. 

Other systematic uncertainties are due to the LAT instrument 
response. The uncertainties in the y-ray selection eflftciency are 
estimated to be 10% at 100 MeV, 5% at 560 MeV, and 20% 
above 10 GeV for the IRFs we used here (lAbdo et al .11201 Oeh . 
The Monte Carlo-based PSF used for this study is known to not 
accurately reproduce in-flight data over the whole energy range 
considered. We verified by means of dedicated simulations that 
this does not significantly affect the determination of the gas 
emissivities considered for the discussion. The energy disper- 
sion is routinely neglected in the likelihood fitting of LAT data 
for limitations in computing power: Monte Carlo simulations in- 
dicate that this approximation causes a bias on the order of 10% 
at 100 MeV decreasing to < 5% above 200 MeV. 

4.2. H I emissivity and CR densities 

The Hi emissivity per hydrogen atom relates to the average 
CR de nsity in each of the regions considered. LAT measure- 
ments (lAbdo et al.l l2009dl) show that the Hi emissivity spec- 
trum in the local ISM is consistent with production via electron 
Bremsstrahlung and nucleon-nucleon interactions by CRs with 
a spectrum consistent with that directly measured in the neigh- 
borhood of the Earth. 

The integrated y-ray emissivity > 100 MeV we measure in 
the Cygnus complex amounts to [2.06+0.1 1 (stat.) ^^^^ (syst.)] x 
10"^^ s"^ sr"^. Figure [TT] shows the Hi emissivity spectrum ob- 
tained for the Cygnus complex and compares it with the expecta- 
tions fo r the local interstellar spectrum estimated in lAbdo et aP 
(l2009dh . The latter is compatible (within 10%) with LAT obser- 
vations at mid latitudes in the third Galactic quadrant in the en- 
ergy range 100 MeV-10 GeV ( Abdo et al. 2009d). The spectrum 
is presented for a uniform spin temperature Ts = 250 K; system- 
atic uncertainties due to the H I opacity correction and to the y- 



ray selection efl&ciency are added in quadrature for display. The 
latter give a non-negligible contribution only at energies larger 
than a few GeV. 

The emissivity of atomic gas in the Cygnus region, averaged 
over ~ 400 pc, is consistent with the local emissivity in the 
100 MeV-100 GeV energy range, except for the deviant point at 
10-30 GeV. This can be explained by the difficulty distinguish- 
ing the H I and CO components at high energies, for which the 
y-ray statistics are limited and the hard cocoon source is brighter 
(Fig. [9]). The emissivity spectrum implies that the CR spectra in 
the relevant energy ranges (~ 1-100 GeV/n for nucleons) are 
similar to those measured in the vicinity of the Earth and in- 
ferred from y-ray observations in the nearby interstellar space 
within 1 kpc. 

The variations in average CR densities along the Local Spur 
between the dense Cygnus complex, two segments in the second 
and third quadrants that exhibit - 2 lower surfa ce densities of 
gas (lAbdo et al.ll2010fl: lAckermann et al.ll201 Ibh . and the mid- 
latit ude diff'use medium with a factor ~ 5 lower surface den- 
sity (lAbdo et al.ll2009dh are constrained to within 10% to 35%. 
This is diflScult to reconcile with the idea of a dynamical cou- 
pling between gas and CR densities (e.g. iBertsch et al.lll993l: 
iHunter et"aDl 1 9971) . They are consistent, on the other hand, with 
the small arm-interarm emissivity contrast estimated from L AT 
data in the third Galactic quadrant (lAckermann et al.|[201 Ibl) . In 
spite of the high column densities of gas, exceeding 10^^ atoms 
cm"^ over many directions within the Cygnus complex, we find 
no hints of excluding CRs from the densest parts of the atomic 
clouds. 

Owing to the bright foreground of the Cygnus complex 
and individual sources, studying the gas emissivity in the outer 
disk of the Milky Way in detail is beyond the scope of this 
study. However, the ratio of the integrated Hi emissivity of 
the outer region over that in the Loca l Spur is (90 + 7)% , in 
very goo d agreement with the results by I Abdo et al.l (1201 Of) and 
lAckermann et al.l (l2011bl) . It confirms in another direction the 
presence of high CR densities beyond the solar circle. 

Located at a distance of ~ 1.4 kpc and / = 80°, the Cygnus 
complex lies at 7? ^ 8.4 kpc from the Galactic center, at a slightly 
smaller radius than the solar system. We measure in this direc- 
tion an emissivity that is consistent with other values found in 
the Local Spur and in three outer segments of the Milky Way 
with Galactocentric radii up to ~ 15 kpc. It yields a decrease 
in H I emissivity < 60% over ~ 6 kpc in Galactocentric radius 
across and beyond the solar circle. Nevertheless, those measure- 
ments span a narrow range in azimuth around the Galactic center 
and may not be representative enough for comparison with ax- 
isymmetric propagation models in order to study the CR gradient 
across the solar circle . 

4.3. CO-bright molecular gas 

If molecular and atomic gas are illuminated by the same CR 
fluxes, we expect the emissivity per hydrogen molecule to be 
twice the emissivity per hydrogen atom, so we can calibrate the 
Xco ratio. We performed a linear fit, qco,i = ^ + 2Xco ' qni,i, 
taking the uncertainties on both emissivities into account to de- 
rive the best linear relation shown in Figure [121 We also give 
the residuals in units of standard deviations. A good linearity is 
found in the 0.1-10 GeV energy range. The highest energy (low- 
est emissivity) points show < 3cr excess of emission associated 
with CO with respect to the Xco ratio determined at low ener- 
gies. The high CO emissivity recorded at 10 - 30 GeV (second 
point) corresponds to a low emissivity in Hi (Fig.[TT]) and may 
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Fig. 12. Top: emissivity per Wqo intensity unit versus emissivity 
per hydrogen atom in the Cygnus complex (for Ts = 250 K). 
The points correspond to the different energy bins; the emissivi- 
ties decrease with increasing energy. The red line gives the best 
linear fit taking uncertainties on both axes into account. Bottom: 
residuals in units of standard deviations as a function of H I emis- 
sivity. 



result from a fluctuation in the difficult spatial separation be- 
tween the atomic and molecular components when photons are 
sparse given the hard 2° source that partially overlaps the CO 
peaks. Up to 10 GeV, the linearity is good, so there is no sign 
of CR exclusion from the dense cores of this giant molecular 
complex. 

The slope of the best-fit linear relation provides a value 
of Xco = (1.68 ± 0.05) X 10^0 cm'^ (K km s-^)-^ in the 
case of Ts = 250 K. We obtain Xco = (1.58 ± 0.04) x 10^^ 
cm"^ (K km s"^)"^ in the limit of small Hi optical depth and 
Xco = (2.55 ± 0.08) X 10^^ cm'^ (K km s-^)-^ in the case of 

= 100 K. The uncertainties in qm associated with the H I spin 
temperature are particularly severe for the high-density clouds 
of the Cygnus complex. High optical depths (low spin temper- 
atures) imply a large increase in N(tl l), therefore substantially 
lower CR densities. Given the y-ray luminosity of the molecular 
clouds, this subsequently implies a significant increase in their 
estimated masseo The systematic errors on the y-ray selection 
efl&ciency cancel out to the first order in the estimate of the Xco 
ratio. 

The conversion factor Xco = [1.68 ± 
0.05 (Stat.) !^:^^(Hi opacity)] x lO^^ cm'^ (K km s-^)-^ 
is consistent with other LAT measurements in the Local Spur, 
which r ange from 1.5 to 2 xlOj^^m""^ (Kkms-i)-i (lAbdo et al.l 
1201 Ofl: lAckermann et al.l l2011bh . From the diff'erent y-ray 
measurements in the Galactic plane, the Xco ratio at the solar 



circle appears well defined. It is, however, significantly greater 
than in nearby well-resolved clouds off' the plane in Cassiopeia 
and Cepheus ( Abdo et al. 20101). Whether the discrepancy is 
due to the sampling resolution or to an intrinsic Xco variation on 
diff'erent scales inside a cloud will be investigated in the future. 

Using the Xco ratio, we estimated the CO-bright molecular 
mass in the complex. For this purpose we considered the region 
at 74° < / < 86°, -5° < b < 8°, where most of the gas as- 
sociated with the Cygnus complex is located. Assuming a dis- 
tance of 1 .4 kpc and a mean atomic weight per hydrogen atom 
in the ISM of 1.36, we obtain a mass 2.3 ^ x lO^MO (where 
the uncertainties are dominated by the H I opacity correction). 
This value (taking the diff'erent assumpt ion on the distance into 
account) is consistent with the results by lSchneider et al.l (l2006l) 
based on higher resolution, multi-isotopolog CO observations, 
and it depicts Cygnus as a super-massive molecular complex. 

The small amount of CO-bright molecular gas in the outer 
region of the Milky Way in this longitude window (Fig.|4]) makes 
determining its emissivities extremely sensitive to the details of 
the model (including point sources), so we do not consider it for 
scientiffc interpretation. 

4.4. Dark neutral gas 

Including the Ay excess map in the model corresponds to an 
increase of 250.6 in the logarithm of the likelihood (for ten ad- 
ditional degrees of freedom). This corresponds to a signiffcant 
detection of y-ray emission associated with Ay excesses, for- 
mally equivalent to a ~ 2lcr conffdence level. Figure [T3l shows 
the emissivity per Av,exc unit, ' versus the emissivity per hy- 
drogen atom, qm, in the Cygnus complex. A good linear correla- 
tion is found between the two emissivities over three decades in 
energy, proving that y-ray emission associated with Ay excesses 
comes from the same physical processes as that associated with 
H I. Ay residuals therefore trace interstellar gas. 

With a procedure analogous to what is adopted to estimate 
^co, we can use the emissivity per hydrogen atom to cali- 
brate the dust-to-gas ratio in the dark neutral phase Xa^ = 
A^(H)/Av,exc. We obtain Xa^ = (28 ± 2) x 10^^ cm'^ mag-^ 
in the case of Ts = 250 K, Xa^ = (48 ± 3) x 10^^ cm'^ 
mag-i in the case of Ts = 100 K and Xa^ = (27 ± 2) x 10^^ 
cm"^ mag"^ in the case of optically thin medium, so Xa^ is 



l20 2 j^^g 



The same level of uncertainty would affect the Xco derivation from 
another total gas tracer such as the dust column-density. 



[28 ± 2 (Stat.) !f (H i opacity)] x 10^ 

Assuming a standa rd total-to- sele ctive extinction ratio Ry = 
Av/E(B - V) = 3.10 (IWegnedl2003l) . the dust-to-gas ratio just 
estimated appears to be ~ 50% higher than the average value 
in the diffuse ISM, A^(H)/E(B - V) = 58 x 10^^ cm'^ mag-^ 
(Bohlin etaLlll978|), and a factor of three higher than what is 
inferred for the dark phase in local clouds from y -ray measure- 
ments MH)/E(B- V) 30 X 10^^ cm-2 mag-i (iGrenier et al.l 
120051: 1 Abdo etal ] |2010j) . The discrepancy is confirmed by the 
extinction data. By fitting the latter with the H i and CO maps 
(§[223, we obtained A^(H I)/ Av = 29.6+0. 1 x 10^^ cm'^ mag-^ 
(statistical error only for Ts = 250 K). A possible explanation 
is provided by an anomalo us Ry ratio driven by a peculiar dis- 
tribution of dust grain sizes (ICardelli et al.lll989h . lStraizys et al.l 
(1999) indeed report an anomalous extinction law in the Cygnus 
region, showing stronger extinction in the violet and near UV 
region. 

The chemical state of the dark neutral gas cannot be deduced 
from y-ray observations. Whereas there are compelling theoret- 
ical and observational reas ons to believe that CO-quiet H2 is 
ubiqu i tous in the ISM (e .g. IWolfire et 10120101: iMagnani et al.l 
I2OO3I: iLanger et"an l2010l) . we cannot exclude that part of the 
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Fig. 13. Top: emissivity per Av,exc unit versus emissivity per hy- 
drogen atom in the Cygnus complex (for Tg = 250 K). The 
points correspond to the different energy bins; the emissivities 
decrease with increasing energy. The green line gives the best 
linear fit taking uncertainties on both axes into account. Bottom: 
residuals in units of standard deviations as a function of H I emis- 
sivity. 



dark gas traced by Ay excesses is missing cold atomic gas, es- 
pecially since the dark neutral phase appears at the interface be- 
tween the atomic a nd CO-bright phases in the nearby clouds 
(iGrenier et al.ll2005"b . Temper atures as low as 40 - 7 K were 
measured in cold H I clouds (iHeiles & Trolandll2003h . and self 
absorption can be strong when cold clouds are seen against more 
diff'use warm H I. The Ay excesses in Fig. [5] partially overlap an 
H I self-absorption feature associated with the Cygnus complex 
(iGibson et aLlhoOSL Fig. Id). Nevertheless, the Hi to H2 transi- 
tion is very dynamical, both in space and time, and it is difl&cult 
at this stage to conclude anything about the exact mix of cold 
dense H l and diff'use CO-quiet H2 that forms the dark neutral 
phase in the outskirts of CO-bright molecular clouds. 

Regardless of its nature, the mass of the dark neutral gas in 
the Cygnus complex at 1.4 kpc amounts to 0.9 x lO^MO. 
Adding an atomic mass of 5 x lO^MO and including the CO- 
bright mass estimated above the total interstellar mass of the 
Cygnus complex amounts to 8 x lO^MO. 

Assuming that all the dark neutral gas is molecular, we 
can calculate the molecular dark-gas fraction /dg = (^moi - 
Mco)/^moi, which amo unts tcB 0.27 ± 0.02 , in excellent agree- 
ment with the model bv lWolfire etaD(l2010b. The dark-sas frac- 
tion is also consistent with the one by lAbdo et al.l (l2010j) for the 
nearby Cepheus and Cassiopeia clouds, which have a factor of 
two lower column densities and masses < 2% of that contained 



The dark neutral gas fraction is very stable against the choice of 
H I spin temperature, therefore the error on the dark-gas fraction is only 
statistical. 



in the Cygnu s comp lex. This also agrees with the prediction by 
IWolfire et al.l(l2010l) that the dark-gas fraction is fairly indepen- 
dent of the mean cloud column density and total mass for giant 
molecular clouds. 



5. Conclusions 

We performed an analysis of y-ray emission across the entire 
Cygnus region measured by the Fermi LAT in the energy range 
100 MeV-100 GeV. We built a general model for the region able 
to satisfactorily reproduce the LAT data. The model includes 
extended sources that have been detected over the interstellar 
emission model described here in association with the Cygnus 
Loop and y Cygni supernova remnants and with a cocoon of 
freshly-accelerated CRs in the innermost part of the Cygnus X 
region. They are discussed in detail in companion papers. 

We measured the average Xco = A/^(H2)/Wco 
factor for clouds in Cygnus, finding a value [1.68 ± 
0.05 (Stat.) !^:^^(Hi opacity)] x lO^^ cm'^ (K km s'^)"! 
that is well consistent with other LAT measuremen ts for cloud 
complexes in the Local and Pe rseus spiral arms (lAbdo et al.l 
l2010j: IXckermann et al.l [201 Ibh . These Xco ratios, averaged 
over complexes, are, however, significantly higher than the 
values found at higher sampling re solution in nearby clouds of 
the Gould Belt (lAbdo et al.ll2010fl) . Thanks to the correlation 
between dust and y-ray emission excesses, we detected the 
presence of conspicuous masses of dark neutral gas not traced 
by the combination of the Hi and CO lines, with total mass 
~ 40% of the mass of the clouds traced by CO. The good 
correlation over three decades in energy between the y-ray 
emissivity per Ay excess unit and per H atom strengthens the 
interpretation of such excesses as produced by dark neutral 
gas. The neutral gas in the Cygnus complex, combining atomic, 
CO-bright, and dark masses, amounts to 8^^ x lO^MO at a 
distance of 1.4 kpc. 

The emissivity of atomic gas measured over the whole 
Cygnus complex is consistent with other estimates in the local 
interstellar space. We do not find evidence of any possible ex- 
clusion of CRs by enhanced magnetic fields in the dense clouds. 
The emissivity per hydrogen atom compares with LAT estimates 
in other regions of the local and outer Galaxy, regardless of dif- 
ferences in gas surface density by about one order of magnitude 
and in Galactocentric radius by ~ 6 kpc. This uniformity does 
not support models based on the dynamical coupling of CRs with 
matter densities or predicting a strong emissivity gradient toward 
the outer Galaxy. 

The CR population averaged over the whole Cygnus com- 
plex (~ 400 pc) is similar to the Local Spur average, in spite 
of the embedded regions of conspicuous massive star-formation 
and potential CR accelerators. Their impact on the CR popu- 
lation is only detected in the innermost region bounded by the 
ionizati on fronts from the massi ve stellar clusters on a scale < 
100 pc (lAckermann et al.|[201 lah . No counterpart to the broadly 
distributed excess of y-ray emission seen at energies > 10 TeV 
at 65° < / < 85° ( Ab do et al...2007.,2008.) have been detected at 
GeV energies so far. 
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